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摘 ， 要 :针对 深海 采矿 系统 长 距离 重 直 输送 管道 顶端 不 同 约束 条 件 对 复杂 流 固 耦 合 振动 特性 的 影 
响 ,采用 基于 双向 流 固 耦合 的 模 态 分 析 方法 ,重点 研究 了 匀 接 约束 和 弹性 约束 对 管道 国有 频率 及 振 
型 的 影响 ;结合 模 态 分 析 结 果 以 及 管道 外 壁 所 受 应 力 和 变形 情况 ,对 弹性 约束 管道 进行 谐 响应 分 
析 , 并 根据 分 析 结 果 对 结构 进行 了 优化 设计 。 研 究 表明 :对 比 两 种 约束 条 件 模 态 分 析 结 果 , 弹 性 约 
) 束 管道 固有 频率 最 大 、 振 幅 最 小 , 且 振 型 最 为 简单 ,实际 工程 中 管道 顶端 选用 弹性 约束 最 合适 ; 当 激 
于 励 频 率 处 于 第 7 阶 固有 频率 附近 时 ,管道 的 振动 响应 最 为 激烈 ;通过 增 焊 强 环 圈 的 方法 对 输送 管道 
DD 进行 优化 ,对 比 优化 前 后 模 态 分 析 结 果 , 证 明了 优化 的 可 行 性 。 
“关键 词 :输送 管道 ;边界 约束 ;固有 频率 ; 模 态 振 型 ; 谐 响应 分 析 ; 结构 优化 
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Abstract :In view of the influence of different constraints on vibration characteristics of transportation pipe- 
line ,the research mainly focuses on the effects of hinge and elastic constraints on the natural frequencies 
and modes of the pipeline by the bidirectional fluid structure interaction. With the results of modal analysis 
and the stress and deformation on the pipe outer wall ,the harmonic response of elastic-constrained pipe is 
analyzed,and the structure is optimized according to the analysis results. The results show that the elastic 
constraint pipe has the largest natural frequency ,the smallest amplitude ,and the simplest vibration mode; 
therefore ,it is the most suitable choice for the top of the pipe in the actual project. When the excitation fre- 
quency is close to the seventh-order natural frequency ,the vibration response of pipeline is the most in- 
tense. The method of welding strong ring is used to optimize the pipeline ,and the results of modal analysis 
before and after the optimization are compared ,which proves the feasibility of the optimization . 
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完整 的 海底 矿产 资源 开发 系统 ,从 水 面 至 海底 ， 
主要 由 大 型 水 面 支持 船 . 深 水 悬臂 式 长 距离 输送 管 
道 , 水 下 大 功率 提升 泵 和 海底 集 矿 装置 等 部 分 组 成 。 
其 中 ,长 度 上 千 米 的 长 距离 垂直 输送 管道 是 整个 系 
统 中 最 为 薄弱 ,技术 难度 最 大 的 部 分 ,也 是 最 为 关键 
的 核心 内 容 。 管 道外 部 承受 海流 波浪 等 环境 荷载 
的 作用 ,内 部 受到 有 一 定 流 速 的 矿 浆 影响 ,在 内 外 流 
联合 作用 下 会 产生 流 固 耦 合 振动 。 当 外 界 激 励 频率 
接近 管道 固有 频率 时 ,会 引起 管道 的 共振 ,影响 管道 
的 工作 稳定 性 和 疲劳 寿命 。 因 此 ,深海 采矿 长 距 
离 输送 管道 的 振动 特性 研究 ,对 结构 的 性 能 分 析 、 优 
化 设计 具有 重要 意义 。 

深海 采矿 长 距离 输送 管道 振动 特性 分 析 方 法 主 
要 包括 理论 分 析 、 数 值 模拟 以 及 试验 研究 ”” 。 其 
中 5S 朱 卫 平等 中 通过 Hamilton 原理 推导 出 了 外 部 横 
疡 猎 励 作用 下 固 - 匀 支 撑 管道 振动 微分 方程 ,并 用 
Calerkin 法 求解 验证 了 振动 方程 的 正确 性 ; 俞 树 荣 
等 四 推导 出 海洋 立 管 耦合 振动 方程 ,并 通过 Hermit 
搬 介 函数 法 对 振动 方程 进行 离散 求解 Heili7 通过 
乱 项 迁 代 法 研究 了 大 位 移 流体 结 者 构 相互 作用 问题 的 
合家 合 ( 单 片 ) 解 ; Huang 等 5 通过 消 元 Galerkin 法 
研究 了 不 同 边界 条 < 件 下 输送 管道 固有 频率 ,得 到 了 
相同 边界 条 件 下 的 输送 管道 固有 频率 方程 ; Dai 
等 “采用 传递 矩阵 法 ,通过 引入 动态 刚度 矩阵 ,计算 
下 简单 直 管 和 弯 管 系统 的 固有 频率 ;文献 [10-11] 研 
内 外 流 因 耦合 作用 对 管道 自由 端 偏 移 量 的 影 
响 S 匡 兵 等 中 将 两 端 固 支 薄 壁 圆柱 管 简 化 为 Timo- 
shEiko 梁 模型 ,自制 试验 装置 研究 薄 壁 圆柱 管 的 共 
振 频 率 和 振 型 ;Kang 等 "研究 了 在 均匀 流 场 环 境 
下 ,张力 对 长 圆柱 体 涡 激 振动 响应 的 影响 ; 文献 
[14-15] 对 比 了 柔性 管道 的 数值 计算 和 物理 模型 试 
验 , 指 出 物理 模型 试验 对 于 理解 管道 复杂 动力 响应 
的 必要 性 。 

以 上 学 者 进行 研究 时 大 多 将 输送 管道 简化 为 其 
臂 梁 ,而 在 真实 工作 环境 中 管道 与 采矿 船 还 有 铵 接 
约束 .弹性 约束 等 连接 方式 ,顶端 不 同 约束 条 件 对 应 
管道 模 态 影响 有 所 不 同 。 本 研究 以 深海 采矿 系统 顶 
端 与 采矿 船 的 连接 方式 为 研究 对 象 ,采用 三 维 势 流 
理论 及 修正 的 Morison 方程 推导 内 外 流 作用 下 ,深海 
采矿 长 距离 输送 管道 流 固 耦合 振动 方程 ;同时 ,利用 
有 限 元 方法 ,计算 顶端 球 贸 接 和 弹性 连接 两 种 不 同 
约束 条 件 下 ,深海 采矿 输送 管道 固有 频率 及 振 型 , 研 
究 工程 中 较 适合 的 连接 方式 ;并 尝试 对 结构 进行 优 
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化 ,中 在 为 实际 工程 中 减 小 管道 振动 对 输送 管道 的 
影响 ,保证 结构 可 靠 性 提供 一 定 参考 依据 。 


1 水 下 输送 管道 的 自由 振动 


1.1 长 距离 输送 固 耦 合 振动 方程 


深海 采矿 系统 简化 模型 如 图 1 所 示 。 管 道 流 固 
耦合 动力 学 计算 简 图 如 图 2 所 示 。 


一 水 面 船 


图 1 深海 采矿 系统 简化 模型 
Fig.1 Simplified model of deep-sea mining system 
Z=0 


y 


llll 


图 2 管道 流 固 耦 合 动力 学 计算 简 图 


Fig.2 Computational diagram of pipe 


fluid-solid coupling dynamics 
假定 计算 流体 为 无 络 、 无 旋 、 有 势 的 理想 流体 ， 
则 对 于 结构 物 所 处 流 场 存在 速度 势 函 数 ,满足 
Laplace 方 程 , 即 
TD ,1 | yoo (1) 
or rr or rr 09090 0z 
式 中 :BD(r,0,z,i) 为 流体 速度 势 函 数 ;r、9 为 管道 截 
面 的 极 坐标 。B(r,0,z,it) 满 足 如 下 边界 条 件 。 
1) 在 固定 界面 上 ,有 


0 (2) 
2) 不 考虑 重力 的 影响 ,在 液体 表面 ,有 
D(r,0,z,t) =0,(z=0) (3) 


3 ) 在 流 场 无 穷 远 边界 处 ,有 
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BD(r,0,z,t) =0,(7 一 oo) (4) 
4) 当 909=0 时 ,有 
OP(r,0,2,t) P(rT,T,2,t) 
00 00 Sy 3) 
在 流 固 耦合 面 上 ,管道 的 振动 速度 与 流体 的 振 
动 速度 相同 ,有 
2 =z) cos0 (6) 
式 中 ,w_,(z,t) 为 管道 来 流 方向 速度 。 
ww,_(2,t) =U(2) 字 (7) 


式 中 U(z) 为 管道 在 海水 中 的 振 型 函数 。 
将 边界 条 件 式 (2) ~ 式 (5) 代 入 式 (1) ,得 
hy J (Air) cos( Aiz) 
L di Aa AiJ' (Airo) 


Dl(r,0,z,t) 一 


fv(8) eos( 8) de (8) 


寺 w 为 管道 在 海水 中 的 固有 频率 。 
P(r,0,z,t) 为 流体 动 水 压力 , 且 
P(r,0,z,t) = -pe (9) 
6 在 流 固 焕 合 面 上 ,有 =r。, 且 
到 2 到 他 万 ( Le 4) 
J UC8) eos( iB) dB (10) 


;只 考虑 管道 在 海水 中 做 自由 振动 ,所 有 变量 都 
按 距 弦 方 式 随时 间 变 化 ,因此 可 以 在 方程 中 将 有 关 
| 中 继 站 外 形 尺 寸 相对 数 千 米 
扬 砍 硬 管 较 小 , 只 考虑 其 质量 对 扬 矿 便 管 的 影 
响 ] ,因此 管道 在 水 中 的 振动 方程 为 

EI Oru(z,t) Tu(zt) 1 
Oz ! 05 
i (st) py Oa) 0 (11) 
05 
将 式 (7) 式 (10) 代 入 式 (11) ,得 
aU(z) (m+ m2 ) w” 
es + EY U(z) + 
2 Ji(A,r)cos(A,z) 
2 A (Airo) 


PE30,z,t) = 


2p,w*cosb 
人 


| v8) cos( hi8)d8 =0 


(12) 


式 中 ;mi 为 管道 质量 ;m, 为 管道 自由 端 集中 质量 ;p， 
2 管道 


为 流体 密度 ;p, 为 管道 材料 密度 ; 27 为 管道 抗 弯 刚 
度 , 令 


m+m,=M (13) 


全 4 不 甘 日 工 || 
| 六 万] 


9 
J UC8) eos( 8) ds (14) 
可 将 式 (12) 写 成 
BISU) ro (M+M* ) Uz) =0 (15) 


式 中 M "为 流体 附加 质量 。 因 此 扬 矿 系 统 输 送 管道 
在 海洋 中 自由 振动 的 湿 模 态 ,可 以 看 作 是 管道 自 里 
及 其 周边 流体 动 水 压力 附加 到 管道 上 产生 的 振动 。 
求解 (12) 式 得 

U(z) = Dia(kz) + D,B, (kz) + Dyy (kz) + 


,GDas (ks) + D,B,(hz) +D, (Yilie) + + a)* 


D, (a.( hz) re) cos(A,z) (16) 


其 中 
i -+ cos( kz) 
ra _ cosh( hz) — cos(kz) 
| 。 . (17) 
yy) = sin( kz) 
as =- sin( hz) 
. cosh(kz) cos(kz) 
Nisin | Br Ra 
Qi( kz) = 7 
cosh( kz) + hz) 
Ajsin( A,L) | 及 | 
Bil( hz) = 7 
(18) 
Ph [+ | 
yi(hz) = 7 
oo) -~ 
Nsin(AL) | -| 
i k CE 
Ce -2piTrow” 
要 4 J (Airo) 
LA EN! —M)w We +p;Trow | 
(19) 
4 _ Mo” 
此 兰 ET (20) 


利用 振动 方程 式 (7)、 式 (14) 和 式 (16 ) ,结合 
表 1 不 同 约束 条 件 下 管道 的 边界 条 件 , 即 可 求 其 固 
有 频率 和 振 型 。 匀 接 约束 的 形式 采用 球 贸 接 。 


表 1 不 同 约束 条 件 下 管道 边界 条 件 


Tab.1 Boundary conditions of pipelines 
under different constraints 
约束 匀 接 约束 弹性 约束 
(0,1) =0, = 
约束 端 92u(0,1) 
El— =0 9 Ou(0,D)] 
EE | |= fuC 0,1) 
2 
BI ou Lt) -0， 二 
自由 端 0z 0 2Z 
1 2 2 
0 ou(L,t) 0 ou(L,i) 
一 | EI—— |= -一 | EE/— oO |= 
| oz | 0 | O22 | 
Aisin(A.L 
hr = 1+ DG k=1+ 6 
i=1 i is1 


hk = a(kL) - Y Gia hL) cos( XL) 


< 求解 式 (22 ) , 令 行 列 式 为 0, 行 列 式 中 只 包含 未 
知 疾 ,因此 可 以 求 出 输送 管道 在 海洋 环境 中 的 固 
有 预 率 。 将 上 代入 式 (21) 即 可 求 出 访 的 值 。 将 
种 态 的 值 代入 式 (16) , 即 可 求 出 输送 管道 在 深海 环 
境 钙 自由 振动 的 振 型 。 


1@ 输送 管道 水 平 液 动力 分 析 


,深海 采矿 长 距离 输送 管道 外 部 所 受 的 液 动力 由 
拖 蝎 力 和 惯性 力 组 成 ,其 中 拖 电力 为 由 水 质点 运动 
速度 引起 的 对 管道 的 作用 力 ,惯性 力 为 水 质点 运动 
加 速度 引起 的 对 管道 的 作用 力 。 由 于 管道 直径 与 人 
射流 的 波长 相 比 尺度 较 小 ,因此 适合 用 Morison 方程 
计算 液 动力 。 假 设 波浪 与 海流 同 向 , 均 为 坐标 轴 x 
轴 正 向 , 则 波 流 联合 作用 在 输送 管道 上 的 水 平 液 动 
力 为 
F, = far = | poCoD,(w +v.) |v,+v.| dz+ 


nD? 
4 PwCuvode ( 23 ) 
0 = 人 ereosg (24) 
VR = sing (25) 
12 
， =0. 15 +0.622x| 全 (26) 


式 中 :pw 为 海水 密度 ; C， 为 拖 点 力 系 数 ; Cy 惯性 力 
系数 ;v 为 波浪 水 质点 水 平 速 度 ;v; 为 波浪 水 质点 


k,, = a(kL) =- > GBicos( AL) 


k= e(kL) +2 GelkL)sin(AL) hy = ykL) +2》 Gy(kL)sin(AL) ks = a(kL) + Giaicos(AL) 
ee— Fa i=1 ja 


1.2 方程 求解 
以 管道 顶端 球 铵 接 为 例 ,对 式 (16) 求 导 并 将 边 
界 条 件 代 入 得 


Dk + Dk + Dskis + Dik =0 
Di + Dk + Dskys +D =0 
Di + Ds,ks, + Dskss +D =0 
D,=0 

其 中 刚度 系数 组 成 的 行列 式 为 


(21) 


Aik sin( AL) 加 a 
Py hs 


ks = e(kL) — py Gyicos( AL) (22) 


水 平 加 速度 ;v, 为 海流 水 质点 的 水 平 速度 。 


2 有 限 元 建 模 及 模 态 分 析 


2.1 输送 管道 的 有 限 元 建 模 


针对 分 布 在 水 深 800 m 以 下 的 富 销 结 核 区 和 富 
锰 结 核 区 ,本 研究 基于 ANSYS Workbench 有 限 元 分 
析 软 件 , 以 计算 总 长 为 1000m 的 长 距离 重 直 输送 管 
道 建立 计算 模型 。 模 型 由 输送 管道 及 中 继 站 组 成 ， 
如 图 3 所 示 。 外 部 流 场 域 与 内 部 流 场 域 均 采用 流体 
单元 (Fluid 30 ) , 输送 管道 采用 结构 单元 (Solid 
186) ,在 管道 的 内 外 壁 定义 无 滑 移 流 固 耦合 边界 。 
为 得 到 精确 的 计算 结果 ,管道 及 内 外 流 场 域 均 采用 
六 面体 网 格 , 外 部 流体 网 格 按照 越 靠近 管道 处 越 密 
的 原则 划分 ,网 格局 部 图 与 网 格 剖 面 图 如 图 4 图 5 
所 示 ,结构 网 格 单元 数量 为 1 656 239 ,流体 网 格 单元 
数量 为 4830753。 计 算 域 模型 如 图 6 所 示 。 


多 3 输送 管道 几何 模型 
Fig.3 Geometric model of the conveyingpipe 
假设 在 作业 过 程 中 ,采矿 船 与 输送 管道 的 运动 
相对 较 小 ,将 管道 顶端 考虑 成 匀 接 约束 或 者 弹性 约 
束 ,由 于 扬 矿 软 管 的 刚度 较 小 , 且 附 着 有 浮力 材料 ， 
对 扬 矿 便 管 的 作用 力 可 以 忽略 不 计 , 可 将 中 继 站 下 
端 作为 自由 端 来 处 理 。 以 4 级 海 况 为 例 ,输送 管道 
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计算 参数 为 :输送 管道 总 长 1000 m, 内 径 420 mm, 外 
径 500 mm, 壁 厚 40 mm ,密度 7 850 kg/mi ,单位 长 度 
质量 355. 1 kg/m, 弹性 模 量 210 GPa; 海水 密度 
1037 kg/m ; 矿 浆 密度 1 250 kg/m ;中 继 站 质量 
10t; 海 浪 浪 高 2.5 m, 周 期 10 s。 


图 4 计算 域 网 格 模型 图 


Fig.4 Computational domain mesh grid model 


图 5 模型 网 格 剖 面 图 
Fig.3 Model grid profile 
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2.2 模 态 分 析 


在 ANSYS Fluent 中 定义 外 流 场 域 与 内 流 场 域 ， 
采用 RNG hk-e 满 流 模型 封闭 计算 ,在 外 流 场 域 中 由 
UDF 导入 沿 x 轴 正 方向 的 波浪 水 质点 水 平 速度 v、 
海流 水 质点 的 水 平 速度 v., 内 部 矿 浆 流速 5 m/s, 计 
算 输 送 管 道内 外 流 场 情 况 。 将 流 场 计算 结果 导入 到 
Workbench 中 的 Static Sturctural 以 及 Model 模块 中 ， 
将 流体 对 管道 的 作用 以 预 应 力 的 方式 施加 到 管道 表 
面 , 研 究 不 同 约束 ( 饼 接 、 弹 性 连接 ) 对 输送 管道 固 
有 频率 及 振 型 的 影响 。 计 算 过 程 中 管道 壁面 采用 无 
滑 移 壁面 ,不 考虑 管道 与 流体 的 热 交换 。 由 于 低 阶 
模 态 在 结构 振动 时 起 主要 作用 ,本 研究 主要 研究 前 
12 阶 模 态 。 通 过 Workbench 中 的 Transient Structure 
模块 进行 瞬 态 结构 动力 学 看 合计 算 , 基于 模 态 分 析 


ChinaXiv 合 作 期 刊 


结果 及 瞬 态 动力 学 分 析 绪 果 进行 谐 响应 分 析 , 作为 
动 强度 、 动 刚度 的 校 核 指 标 。 


3 “分析 与 讨论 


3.1 顶端 不 同 约束 条 件 对 管道 模 态 的 影响 


深海 长 距离 输送 管道 顶端 与 采矿 船 有 多 种 约束 
方式 ,本 研究 重点 研究 在 外 部 波 流 联合 作用 以 及 内 
流 作 用 下 , 铵 接 约束 和 弹性 约束 两 种 约束 条 件 对 输 
送 管道 固有 频率 及 振 型 的 影响 。 匀 接 约束 采用 顶端 
球 铵 接 的 约束 形式 ,弹性 约束 通过 在 扬 矿 管 项 端 设 
计 纵 向 的 弹簧 减 振 器 施加 。 根 据 管 道 振动 微分 方 
程 生 ” , 弹 得 减 振 器 的 刚度 设置 为 9 x 10" N/m。 取 
内 流 流速 为 5 m/s, 内 流 压 强 为 2 MPat”|, 模 态 分 析 
得 到 的 各 阶 固有 频率 如 图 7、 图 8 所 示 。 


0.7 


一 上 一 弹性 连接 约束 


模 态 阶 数 
图 7 不 同 约 束 条 件 下 的 输送 管道 回 有 频率 


Fig.7 Natural frequency of pipeline under 


different constraint conditions 
0.8 


一 和 一 3 自由 度 
一 9 一 4 自由 度 
一 4 一 5 自由 度 


固有 频率 /Hz 
三 


0 2 4 6 8 10 12 


模 态 阶 数 
图 8 不 同 自由 度 的 输送 管道 固有 频率 


Fig.8 Natural frequencies of conveying pipelines 


with different degrees of freedom 
从 图 7 中 可 以 看 出 ,对 于 管道 前 12 阶 固有 频 
率 ,管道 顶端 两 种 约束 方式 对 应 的 固有 频率 变化 趋 
势 一 致 ;管道 顶端 匀 接 约束 时 ,各 阶 固有 频率 均 小 于 


顶端 弹性 连接 约束 时 的 各 阶 固有 频率 ,表明 顶端 固 
接 约束 管道 更 容易 被 外 界 激励 产生 共振 ,这 种 现象 
在 高 阶 更 加 明显 。 弹 性 约束 第 1 阶 固有 频率 超过 饮 
接 约束 第 1 阶 固有 频率 6.1% ,而 弹性 约束 第 12 阶 
固有 频率 超过 匀 接 约束 第 12 阶 固有 频率 10. 8% 。 
这 说 明 随 着 模 态 阶 数 的 增加 , 顶端 约束 条 件 对 管道 
固有 频率 影响 越 大 。 图 8 为 输送 管道 顶端 饺 接 约束 
时 改变 自由 度 对 固有 频率 的 影响 。 通 过 对 比 可 以 看 
出 ,相同 阶 数 下 管道 固有 频率 随 着 自由 度 的 增加 而 
减 小 。 这 可 能 是 因为 系统 刚度 随 着 顶端 自由 度 的 增 
加 而 增 大 ,固有 频率 减 小 。 

由 于 重 根 模 态 对 应 的 振 型 形态 相似 , 仅 相 位 不 
同 , 因 此 取 不 同 顶端 约束 条 件 下 输送 管道 第 1 .3 .5、 
7 琶 摸 态 振 型 ,如 图 9 ~ 图 12 所 示 。 从 图 中 可 以 看 
出 0 两 种 不 同 约束 下 同 阶 次 振 型 形态 基本 一 致 ,但 咯 

同 。 顶 端 饺 接 约束 与 弹性 连接 约束 管道 在 靠近 
约束 处 波峰 幅 值 最 小 , 随 着 管 长 的 增加 ,波峰 处 幅 值 
和 忱 加 增 大 , 旦 顶端 饺 接 约束 管道 增长 趋势 大 于 顶端 
弹 隐 连接 约束 状态 。 
CD 同 种 约束 条 件 下 ,管道 振 型 波峰 处 幅 值 受到 海 
大 洪 度 影响 :在 水 深 0 ~700 m, 振 型 波峰 处 幅 什 随 模 
态 阶 数 的 增加 而 增加 ;超过 700 m 时 , 振 型 波峰 处 中 
值 随 模 态 阶 数 的 增加 而 减 小 。 因 此 ,进行 深度 超过 
700a 的 水 下 作业 时 要 注重 考虑 低 阶 模 态 对 管道 振动 
的 野 响 。 管 道 顶端 无 论 采 用 何 种 约束 ,最 大 变形 均 发 
生 奉 靠近 中 继 站 一 侧 。 对 比 两 种 不 同 约束 条 件 ,顶端 
泗 呈 约束 时 管道 固有 频率 较 高 , 灾 曲 振动 波峰 处 由 
值 较 小 ,实际 工程 中 选用 顶端 弹性 约束 最 合适 。 


9 一 一 1 阶 振 型 
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Fig.9 Mode shapes of hinged connection at the top of pipeline 
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图 10 管道 顶端 弹性 约束 模 态 振 型 


Fig. 10 Mode shapes of elastic connection at the top of pipeline 
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图 11 不 同 顶端 约束 模 态 振 型 对 比 


Fig.11 Comparison of modes with different top constraints 
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图 12 不 同 约束 条 件 下 管道 弯曲 处 波峰 幅 值 


Fig.12 Peak amplitude under different constraints 


波 流 冲 击 下 管道 流 固 耦合 瞬 态 应 力 分 析 
海洋 环境 中 波浪 与 海流 联合 冲击 管道 ,外 部 流 


S12 


场 域 流体 运动 状态 被 改变 ,流体 与 结构 之 间 产 生 泊 
松 耦 合 与 结合 部 耦合 ,对 管道 外 壁 产 生 较 大 作用 力 ， 
导致 管道 发 生 显 车 变形。 内 外流 场 的 边界 条 件 设 置 
为 :内 流 由 下 向 上 运动 ,进口 处 采用 速度 进口 边界 条 
件 , 出 口 处 采用 开放 性 边界 条 件 , 荷载 加 载 于 内 局 
面 。 外 部 流体 的 运动 方向 平行 于 管道 的 径 向 方向 ， 
进口 处 采用 速度 进口 边界 条 件 ,出 口 处 采用 开放 性 
边界 条 件 。 其 他 面 均 采用 对 称 边界 条 件 。 

图 13 ~ 图 15 为 顶端 不 同 约束 管道 在 受到 瞬 态 流 
场 作用 下 的 总 变形 和 等 效应 力 。 从 图 中 可 以 看 出 , 顶 
端 弹性 约束 与 顶端 饺 接 约 束 相 比 , 管 壁 最 大 位 移 减 小 
-0.098 3 ~1.716 m ,等 效应 力 减 小 0.2 ~565 MPa。 


0 0.5 1.0 15 2.0 
时 间 /s 


(a) 最 大 位 移 
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0 0.5 1.0 和 2.0 
时 间 /s 
(b) 最 大 等 效应 力 
图 13 不 同 约束 条 件 下 的 最 大 位 移 和 最 大 等 效应 力 


Fig.13 ” Maximum displacement and maximum 


equivalent stress curves under different constraints 
对 比 两 种 不 同 约束 条 件 ,管道 最 大 位 移 处 与 应 
力 集中 区 域 存 差 异 不 大 。 顶 端 匀 接 约束 管道 与 顶端 
弹性 约束 管道 最 大 应 力 出 现在 管道 最 大 位 移 处 ,应 
力 变化 趋势 基本 一 致 ,但 管道 顶端 为 弹性 约束 时 ,其 
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及 位 移 变 化 较 小 ;顶端 弹性 约束 时 管道 最 大 位 移 及 
最 大 应 力 最 小 ,在 实际 工程 中 应 当 适 当 考 虑 管道 顶 
端 弹 性 约束 。 


ET 


ural- 淆 接 


(a) 总 变形 


(b) 等 效应 力 


图 14 ” 铵 接 约束 总 变形 和 等 效应 力 


Fig. 14 Total deformation and equivalent 


stress of hinged constraints 


(a) 总 变形 


(b) 等 效应 力 


图 15 弹性 约束 总 变形 和 等 效应 力 


Fig.15 Total deformation and equivalent 


stress of elastic constraint 
3.3 弹性 约束 管道 的 谐 响 应 分 析 


基于 深海 采矿 系统 长 距离 垂直 输送 管道 湿 模 态 
分 析 绪 果 与 力学 分 析 结 果 , 采 用 模 态 三 加 法 研究 顶 
端 弹性 约束 管道 在 不 同 持续 频率 的 周期 荷载 作用 下 
的 动力 响应 ,得 到 谐 响应 分 析 结 末 作 为 动 强度 、 动 刚 
度 的 校 核 指标 。 由 应 力 分 析 计 算 可 得 激 振幅 值 为 
0.5 GPa。 根 据 项 端 弹 性 约束 管道 前 12 阶 固有 频率 
值 ,定义 激励 频率 范围 为 0.1 ~1.0 Hz。 按 定义 的 频 
率 范围 等 分 100 份 ,选取 管道 外 壁 与 外 流 接触 的 流 
回 耦 合 面 作 为 重点 监测 位 置 ,分 析 长 距离 垂直 输送 
管道 的 振幅 、 加 速度 、 应 力 、 应 变 与 外 界 激励 频率 之 


最 大 应 力 值 与 最 大 位 移 值 远 小 于 铵 接 约束 。 由 此 可 
以 得 出 ,在 外 部 波 流 联合 力作 用 下 ,管道 应 力 变 化 及 
位 移 变化 主要 集中 在 管道 上 部 200 mm ,其 余 位 置 应 力 


间 的 响应 关系 。 
针对 管道 流 固 耦合 的 谐 啊 应 分 析 可 得 到 振动 幅 
值 \ 应 变 和 应 力 的 频率 响应 ,结果 如 图 16 所 示 。 
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振幅 /mm 


0 0.2 0.4 0.6 0.8 1.0 
固有 频率 /Hz 
(a) 振幅 -频率 响应 


0 0.2 0.4 0.6 0.8 1.0 
固有 频率 /Hz 
(b) 加 速度 -频率 响应 


0 0.2 0.4 0.6 0.8 1.0 
固有 频率 /Hz 
(e) 应力 -频率 响应 


0 0.2 0.4 0.6 0.8 1.0 
固有 频率 /Hz 
(d) 应 变 -者 率 史 应 


图 16 幅 值 ,应力 和 应 变 的 频率 响应 曲线 


Fig. 16 Frequency response curves of amplitude, 


stress and strain 
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由 图 16 可 以 看 出 , 随 着 外 界 激励 频率 的 增 大 ， 
管道 最 大 变形 处 的 振幅 、 加 速度 、 应 力 、 应 变 均 在 第 
1 阶 (0. 156 Hz) ,3 阶 (0.264 Hz) 5 阶 (0.363 Hz)、 
7 阶 (0.459 Hz) .9 阶 (0.555 Hz) .11 阶 (0.649 Hz) 
回 有 频率 附近 出 现 峰 值 。 其 中 振幅 、 加 速度 、 应 力 、 
应 变 的 最 大 峰值 出 现在 第 7 阶 固有 频率 (0. 459 Hz) 
附近 ,第 11 阶 固有 频率 (0. 649 Hz) 附近 出 现 第 2 大 
峰值 ,其 他 模 态 阶 数 附近 响应 峰值 较 小 ,对 结构 影响 
较 微 弱 。 因 此 ,下 文 只 对 当 激励 荷载 为 0.459 Hz 及 
0. 649 Hz 时 管道 响应 进行 分 析 。 对 比 模 态 分 析 结 
与 谐 响应 曲线 ,响应 峰值 频率 与 对 应 固有 频率 的 频 
率 差 第 1.3. 35. 7 .9 、11 阶 分 别 为 1.35% 、0.6%、 
0. 52% .0. 15% .0.74% .0.28% 。 

从 图 16(a) 中 可 以 看 出 :由 于 存在 沿 x 轴 方 向 
的 波浪 和 海流 冲击 ,管道 沿 x 轴 方 向 的 位 移 响应 最 
大 ;在 最 大 峰值 处 ,x 轴 方 向 与 y 轴 方 向 位 移 峰值 为 
54.17 mm 和 8.09 mm, 而 由 于 管道 顶端 被 约束 ,所 
以 在 沿 z 轴 方 向 位 移 响 应 几乎 为 0。 从 图 16(b) 中 
可 以 看 出 :由 于 管道 项 端的 约束 ,x 轴 方 向 加 速度 响 
应 峰值 远 远 大 于 y 轴 及 = 轴 加 速度 响应 峰值 ; 当 激 
励 频率 为 0.46 Hz 时 ,x 轴 方 向 加 速度 响应 峰值 为 
0. 453 m/s ; 当 激 励 频率 为 0. 649 Hz 时 ,x 轴 方 向 加 
速度 响应 峰值 为 0.23 m/s 。 由 此 得 出 结论 ,长度 为 
1 000 m 的 深海 采矿 输送 管道 受到 外 部 激励 时 ,沿海 
流 与 波浪 方向 会 产生 较 大 位 移 和 加 速度 响应 。 在 进 
行 结构 系统 优化 时 可 以 通过 添加 阻尼 器 来 避免 一 定 
范围 内 的 共振 , 减 小 振动 产生 的 变形 。 

从 图 16(c) 和 (d) 中 可 以 看 出 :由 于 存在 慌 直 z 
轴 方 向 上 的 波 流 联合 冲击 ,监测 位 置 z 方向 的 应 力 
应 变 响应 远 远 大 于 其 他 两 个 方向 ; 当 激 励 频 率 为 
0. 46 Hz 时 ,z 轴 方 向 上 的 应 力 响 应 为 18. 16 MPa , 应 
变 响应 为 9.3 x 10 mm; 当 激 励 频 率 为 0. 649 Hz 
时 ,z 轴 方 向 上 的 应 力 响 应 为 5.71 MPa, 应 变 响 应 为 
2.91 x10“”mm。 可 以 通过 安装 振动 传感器 对 管道 
进行 压力 与 速度 的 振动 数据 监测 。 这 些 分 析 可 以 为 
管道 结构 的 设计 与 优化 提供 技术 支持 与 理论 依据 。 


3.4 优化 设计 


通过 模 态 分 析 和 谐 响应 分 析 结 果 可 以 得 知 , 深 
海 采 矿 长 距离 垂直 输送 管道 在 第 7 阶 模 态 振动 响应 
最 为 强烈 , 振 型 以 弯 振 为 主 , 发 生 多 次 弯曲 变形 。 反 
复 大 变形 振动 会 使 管道 产生 疲劳 破坏 ,影响 深海 采 
矿 系统 工作 稳定 性 。 因 此 ,在 深海 采矿 系统 整体 结 
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构 不 变 的 情况 下 ,通过 在 输送 管道 外 壁 增 焊 强 环 圈 管 i 


来 降低 管道 国有 频率 , 避 开 较为 强烈 的 共振 , 如 图 
17 所 示 。 根 据 谐 响应 分 析 可 得 共振 时 结构 变形 最 
大 部 位 相隔 30 mm ,因此 每 间隔 30 m 焊接 强 环 圈 。 强 
环 圈 高 2m, 厚 0.2m。 图 18 ~ 图 20 为 输送 管道 优 
化 后 模 态 分 析 结 


图 17 输送 管道 优化 后 的 部 分 模型 
Fig.17 Partial model of optimized pipeline 
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图 18 输送 管道 优化 前 后 固有 频率 对 比 


Fig.18 Natural frequency comparison before and 
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after pipeline optimization 
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图 19 ”输送 管道 优化 前 后 各 阶 模 态 振 型 最 大 幅 值 对 比 


Fig. 19 Comparison of maximum amplitude of each 


mode shape before and after pipeline optimization 
从 图 18 中 可 以 看 出 :输送 管道 优化 后 各 阶 固有 
频率 较 优 化 前 降低 35. 85% ~ 54. 3% ,有 上 且 各 阶 固有 
频率 之 间 差 值 减 小 ,第 12 阶 固有 频率 为 0.419 Hz， 
避 开 了 响应 最 激烈 的 0.46 Hz 和 0.649 Hz。 从 图 19 


中 可 以 看 出 ,优化 后 管道 各 阶 固有 频率 对 应 振幅 均 
有 所 减 小 ,其 中 最 大 振幅 仅 为 2. 653 6 mm ,为 优化 前 
的 67% 。 从 图 20 中 可 以 看 出 ,优化 后 响应 位 移 大 幅 
降低 ,说明 此 改进 方法 确实 对 管道 抗 振 性 能 起 到 一 
定 作 用 。 


振幅 /mm 
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频率 /Hz 
图 20 输送 管道 优化 后 位 移 -频率 响应 


Fig. 20 Displacement frequency response of optimized pipeline 


4 结 论 


为 了 分 析 顶 端 不 同 约束 条 件 对 深海 采矿 系统 长 
距离 垂直 输送 管道 在 海洋 环境 中 的 振动 特性 影响 ， 
本 研究 采用 理论 推导 和 模拟 仿真 相 结 合 的 方式 进行 
研究 。 首 先 采 用 三 维 势 流 理论 推导 内 外 流 作用 下 深 
海 采矿 长 距离 输送 管道 流 固 耦合 振动 方程 ,后 续 设 
计 了 一 系列 计算 工 况 , 分 析 了 管道 顶端 匀 接 和 弹性 
连接 两 种 约束 条 件 下 输送 管道 固有 频率 和 振 型 的 变 
化 规律 。 通 过 以 上 结果 分 析 得 出 以 下 结论 。 

1 ) 顶端 约束 条 件 对 管道 模 态 影响 较 大 。 管 道 项 
端 匀 接 与 弹性 连接 同 阶 回 有 频率 值 相近 , 振 型 形态 
简单 , 均 以 弯 振 为 主 , 顶端 弹性 连接 振 型 幅 值 小 于 锭 
接管 道 。 因 此 通过 对 比 , 实 际 工 程 中 管道 顶端 选用 
弹性 连接 最 合适 。 

2) 通 过 对 弹性 连接 管道 进行 谐 响 应 分 析 发 现 ， 
由 于 内 外 流 联合 作用 ,激励 频率 在 第 7 阶 附 近 
(0.46 Hz) 和 第 11 阶 附 近 (0. 649 Hz) 时 ,管道 振动 
最 为 激烈 。 根 据 特征 点 的 位 移 加 速度 及 应 力 应 变 响 
应 分 布 规律 ,可 以 看 到 在 垂直 于 输送 管道 流 固 耦 合 
面 x 方 向 上 的 振动 位 移 加 速度 响应 最 为 剧烈 ,输送 
管道 z 方 向 上 应 力 应 变 响应 最 为 剧烈 。 

3) 根 据 模 态 分 析 结 果 及 谐 啊 应 分 析 结 果 , 可 以 
通过 对 共振 时 变形 较 大 部 位 增 焊 强 环 圈 的 方法 来 避 
免 或 者 减 小 变形 ;结构 优化 后 输送 管道 固有 频率 降 
低 了 35.85% ~54.3% , 避 开 了 响应 最 为 激烈 的 外 界 


激励 荷载 ,管道 振幅 仅 为 优化 前 的 67% 。 这 证 明了 
优化 方法 的 可 行 性 ,可 以 在 实际 工程 中 起 到 一 定 参 
考 作用 。 
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